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Abstract 

Black phosphorus is a monoatomic semiconducting layered material that degrades exother¬ 
mically in the presence of light and ambient contaminants. Its degradation dynamics remain 
largely unknown. Even before degradation, local-probe studies indicate non-negligible local 
curvature -through a non-constant height distribution- due to the unavoidable presence of in¬ 
trinsic defects. We establish that these intrinsic defects are photo-oxidation sites because they 
lower the chemisorption barrier of ideal black phosphorus (> 10 eV and out of visible-range 
light excitations) right into the visible and ultra-violet range (1.6 to 6.8 eV), thus enabling 
photo-induced oxidation and dissociation of oxygen dimers. A full characterization of the 
material’s shape and of its electronic properties at the early stages of the oxidation process is 
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presented as well. This study thus provides fundamental insights into the degradation dynam¬ 
ics of this novel layered material. 


Introduction 

The understanding and prevention of oxidation processes underpins many successful technologies 
such as the galvanization of metals, steel manufacturing, and the controlled oxidation of silicon 
wafers, to mention a few examples. Furthermore, the exposure to ambient illumination is behind 
the oxidation and subsequent degradation (also known as weathering) of commercial polymers. 
Two-dimensional atomic materials^ 1 2 like graphene, 3 4 hexagonal boron nitride, 5 6 and transition- 
metal dichalcogcnidcs 7 13 are chemically stable when exposed to light under standard ambient 
conditions. On the other hand, black phosphorus (BP) 14 19 has a remarkable puckered structure 
that anticipates its unique chemical properties. 20 24 Uncapped BP transistors exposed to light and 
to ambient conditions break down after a few hours even though the stack contains a large num¬ 
ber of monolayers, thus implying a complete material degradation that originates from random 
locations at the exposed layer. Capping by an A10 a overlayer 23 25 or by hexagonal boron nitride 
(hBN) 26 29 helps to prevent this acute degradation process. Understanding the oxidation of BP is 
a pressing and relevant problem with deep consequences for the science and engineering of this 
layered material and we establish that, among other mechanisms, intrinsic defects can initiate the 
light-induced degradation of BP: one must lower the reactivity of structural defects to delay its 
photo-oxidation. 

A non-constant height profile is apparent in BP even prior to degradation, 22 23 26 implying the 
existence of intrinsic structural defects 30 33 and a non-zero local curvature through the surface. 
The degradation of uncapped samples leads to a visible increase of these height variations, and 
hence to an increased local curvature 20 23 that manifests as pits, bubbles, and bulges. The prop¬ 
erties of membrane-like materials can be tuned by changes in their local shape, 34 45 and it will be 
shown that the propensity of BP to degrade is ultimately linked to its local geometry. 
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Chemisorption barriers for oxygen dimers (O 2 ) on ideal BP are larger than ten electron-Volts 
(two-hundred and fifty kcal/mol, or a thousand kJ/mol) and cannot be accessed through optical 
excitations in the visible spectrum (390 to 780 nm, or 1.6 to 3.2 eV). We realize, nevertheless, that 
these barriers are largely reduced at intrinsic defects and take on typical values for photoinduced 
chemical reactions. 

Indeed, photoexcitations within 3.7 to 5.0 eV lead to the photochemistry of adenine and aminop- 
urine, 46 setting an energy scale for reactions that could be activated by light. Light introduces suf¬ 
ficient energy to break or reorganize most covalent bonds and enables reactions that are otherwise 
thermodynamically forbidden, given that activation barriers of the order of a few eV (>10,000 
Kelvin) are overcome. 

This article represents a departure from other theoretical works 47 52 that study the degrada¬ 
tion of BP as we focus on the role played by intrinsic defects on this material’s chemistry. This 
mechanistic study begins with a structural and electronic ground state consisting of a planar or 
conical BP monolayer and non-interacting CL molecules. We create excited atomistic (hence elec¬ 
tronic) structures by bringing CL molecules through the BP monolayer to gain insight into the 
energetics involved. We then pursue new equilibrium oxidized structures with ab initio atomistic 
optimizations in which the CL molecules have been chemisorbed onto the BP monolayer, and pro¬ 
vide a thorough geometrical and electronic characterization of these systems towards the end of the 
manuscript. Although our present focus is on the oxidation by CL for a self-contained discussion, 
similar studies can be performed to understand the degradation of BP from water molecules and 
other ambient contaminants. An important component of this study is a tool to analyze the shape 
of two-dimensional materials that we have developed. 32 53 56 


3 




Results and discussion 


The energy for O2 to attach to black phosphorus depends on the local shape 

We first investigate how much energy is needed for O 2 to pierce through a BP monolayer. Two 
finite (aperiodic) structures are considered for this purpose: a planar one containing about 600 
atoms, and a conical one containing about 500 atoms 32 (see Methods). 

The local atomistic geometry is established from five local quantities that indicate how the 
material elongates/compresses (the trace and determinant of the metric Tr(g) and Det(g), where g 
is the metric tensor) and/or curves (through the mean curvature H, and the gaussian curvature K) 
along two orthogonal directions, 32 55 56 and how the local thickness T changes with respect to its 
value To for an ideal crystalline BP monolayer. BP is a buckled material, and To is the distance 
among the two sublayers sq and 52 that make up a monolayer. Extensive details of the discrete 
geometry employed to characterize the shape of BP are given in the Methods section. 32 

Activation barriers E a {z) for the oxidation of BP that take full consideration of the spin- 
polarization of the oxygen dimers are estimated by bringing a vertically-oriented O 2 molecule 
-with its center of mass located at a height z along the local normal- into proximity of the BP 
structure, and crossing BP at the largest possible distance among O and P atoms (see reds dots on 
the structural insets on Figure[T|). These activation barriers are computed with ab initio calculations 
without performing any structural optimization (see Methods). 

As seen in Figure [lji, an oxygen dimer requires about 12 eV to rest at the middle of an ideal 
BP monolayer that has the uniform geometry shown in the upper inset ( Tr(g ) = 1, Det(g) — 1, 
H — 0, K = 0, and t/to = 1) in Figure [jji. The energy barrier is symmetric around z = 0, where 
the dimer’s center of mass coincides with the center of mass of the two buckled BP sublayers]- 32 
Such a high energy barrier makes the chemisorption of O 2 onto ideal BP rather unlikely. 

However, the maximum magnitude of E a decreases significantly within intrinsic defects: for 
instance, the line defect with Burgers vector (0,1) 34 induces a dislocation line containing pen¬ 
tagon/heptagon pairs. 30 The energy range for visible light is within 1.6 and 3.2 eV, corresponding 
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Figure 1: The energy barrier E a for CF to oxidize BP depends on the local shape, and the smallest 
barriers on (a) planar or (b) conical structures occur at intrinsic defects. Atomistic structures and 
the local geometry at one sublayer (s i) are shown as insets. The darker/lighter color seen on Tr(g) 
and Det(g) implies larger/smaller interatomic distances at these defects. The low activation barriers 
facilitated by intrinsic defects are within the visible electromagnetic spectrum, which is highlighted 
by the yellow rectangles within the 1.6-3.2 eV energy range. 


to wavelengths within 390 to 780 nm (yellow rectangles in Figure |T[). The maximum value taken 
by E a as an O 2 molecule pierces BP through the geometrical center of the heptagon (red dot shown 
in the leftmost inset, Figure [jj) is slightly smaller than 6.0 eV and accessible via light-induced 
excitations in the near ultraviolet range. In computing these energy barriers with spin-polarized 
calculations in which the total spin is left unrestricted, the spin of O 2 transitions from a triplet state 
at |z| > 2 A of the center of mass of the BP monolayer, to a singlet state for z < 2 A. This change 
of spin configuration of the optimal structures is in agreement with previous reports. 47 

The reduction of the energy barrier at this intrinsic defect is naturally related to the larger metric 
invariants (Tr(g) and Det(g) seen in the inset of Figure la) that lower the electronic repulsion 
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significantly. Crucially, this is only one of many possible structural defects on BP, as this two- 
dimensional material is known to be polymorphic. 31 33 57 61 

The effect of curvature on E a is studied in a conical structure that contains about 500 atoms and 
acquires its largest curvature and compressive strain at its apex, as indicated by the white tones on 
Tr(g) and Det(g) and the large values of H and K. 32 A dimer piercing the apex (upper-left inset on 
Figure [JJt) encounters a maximum value of E a close to 18 eV, a value further away from the reach 
of optically-induced chemistry. The increase of the barrier arises from the structural compression 
that accompanies the creation of curvature at the apex. 32 One can tell “up” from “down” on a 
conical structure, and this distinction makes E a (z) asymmetric (solid, asymmetric black line on 
Figure [jj>). 

The conical structure employed in Figure [IJ) has a dislocation/disclination axis.’ 32 The charac¬ 
teristic structural reconstruction at the edges of BP 31 57 is a manifestation of the polymorphism of 
this material,^" 61 and we found a metastable local structure along the disclination line of the con¬ 
ical structure with two atoms pulled away from a common bond during the structural optimization. 
We envision that similar defects may originate on BP during growth. This metastable structure 
does not localize electronic states within the semiconducting gap either, 31 a fact that will later 
be discussed in greater detail. The barrier E a is very low at this defect, largely overlapping with 
the energy of photons in the visible range, as depicted by the yellow rectangle on Figure |T]b. An 
oxygen dimer attaching to BP at this intrinsic defect requires an activation energy smaller than 4.0 
eV, which is available through an electronic photo-excitation induced by violet light. The accessi¬ 
ble energy barrier for photo-oxidation at this second defect provides conclusive validation of the 
hypothesis that the photo-oxidation and degradation of BP originates at intrinsic defects. 

Chemisorption and dissociation of oxygen dimers 

Once an O 2 molecule is placed within an intrinsic defect aided by a suitable optical excitation, it 
will initially be expelled from BP with a force that has a normal component proportional to the 
slope of the curves in Figure [I] F z — —dE a /dz , and forces felt by individual oxygen atoms lead 
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Figure 2: Unveiling the dissociation of O 2 at intrinsic defects. The first frames to the left show the 
initial vertical placement of O 2 within the intrinsic defects. Oxygen dissociation at these intrinsic 
defects occurs at the frames highlighted in yellow. Frames at the far right show the optimized 
structures, at which distances among O atoms are 2.96 and 3.08 A for at structure (a) and (b), 
respectively. For reference, the equilibrium distance among O atoms in O 2 is 1.25 A. Full movies 
of the dissociation process are available as Supporting Information. 


to the dissociation of O 2 near BP at intrinsic defects. We demonstrate this dissociation mechanism 
next. 

We registered the atomistic dynamics of O 2 placed initially at rest into the two intrinsic defects 
displaying the smallest energy barriers in Figure [Tjand assuming that the atomistic reconfiguration 
is adiabatic. (Under non-adiabatic conditions, the dimer may collect sufficient kinetic energy to 
cross onto a subsequent monolayer.) Full movies showing this process (movies SI and S2) can 
be found as Supporting Information. The fundamental finding is that the distance among O atoms 
increases from 1.25 up to 2.96 A at the end of the structural optimization in Figure 2a, and up to 
3.08 A at the end of the optimization depicted in Figure 2b: O 2 dissociates at intrinsic defects. 
Additional details follow. 

As indicated previously, the dimer was placed vertically at the onset of the optimization, but it 
is expelled from the defect as the O atoms end up dissociated and lining up horizontally. Distances 
from O atoms to the closest P atoms are equal to 1.53 A and the P-P-0 angles subtended to the 
three closest P atoms are equal to 111°, 116°, and 124°. 

The oxidation mediated by the second intrinsic defect provides further evidence for the fact 
that defects help dissociate O 2 due to a non-symmetric force that pulls O atoms apart. The first 
dissociated O atom in Figure [ 2 }) binds to a single P atom with a bond distance of 1.56 A, quite 
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Figure 3: Energetics at the onset of oxidation. The structures gain between 3.1 and 6.8 eV upon 
absorption of individual oxygen dimers, and release from 1.9 to about 7.0 eV upon oxidation, 
indicating how favorable oxidation is once the initial absorption barriers are accessed with the aid 
of optically-induced excitations. Green diamonds indicate the energies at the snapshots displayed 
in Figure |2j 


similar to the distance seen for the O-P single bonds in Figure [2^. The second O atom binds to two 
P atoms with distances of 1.73 A. We address the energetics as BP oxidizes next. 


The energetics of oxidation 

We now study the successive oxidation of two dimers at nearby structural defects. The two dimers 
are placed far above BP prior to oxidation in order to set a suitable reference energy. As the 
first dimer is placed into the pentagon defect in Figure [3jt, it binds to BP and the system gains 
4.3 eV. The system is allowed to relieve forces (see Methods) and it releases an energy of about 
7 eV in the process, making for a highly exothermic reaction. 48 50 Such a large energy release 
could unleash a continued oxidation process -even in the absence of external illumination- if not 
promptly dissipated. 

Continuing the discussion of Figure [3^, a second dimer is placed into BP at a second heptag- 
onal defect, gaining 3.1 eV as it binds to BP. It is energetically unfavorable for oxygen dimer to 
oxidize consecutive heptagonal defects on this dislocation line, as the relaxed structure containing 
two dimers has an energy larger than the structure in which the second dimer is yet to be absorbed. 
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Oxygens on this second dimer did not dissociate, but the distance between the two P atoms in¬ 
creases from 2.3 A to 3.2 A as a common P-P bond is broken; these two P atoms bind to a single 
O atom with bond distances of 1.87 A and a P-O-P angle of 117°. Oxygens on the second dimer 
remain bounded at a larger separation of 1.48 A. There are no significant rearrangements of the 
first two O atoms upon absorption of the second dimer. 

The process is repeated for two O 2 molecules that oxidize the structural defect that was shown 
in Figures [1J> and [2ji. As seen on Figure [3j?, the energy cost upon absorption of the first dimer is 
3.5 eV, which could be accessible if the system adsorbs light with a wavelength of 354 nm. An 
energy release of 7.0 eV is recorded again as the structural forces are relieved. The second dimer 
seen in Figure^ is unable to dissociate and has a distance of 1.49 A among O atoms, and of 
1.80 A among O-P bonds. The structure gains 6.8 eV when the second dimer is absorbed, which 
may be enabled by an optical excitation in the ultraviolet range, and releases about 7.0 eV through 
the atomistic optimization process to make for a slightly exothermic reaction (see Methods). The 
results from Figure 3 indicate that larger energy gains can be attained when the second dimer is 
not in close proximity to the first absorbed Cb molecule. 

The substantial energy gains upon chemisorption recorded in Figure [3] imply that the oxidized 
structures are more stable than BP containing intrinsic defects, providing a clear picture of how 
reactive the material can be at these defects. Even though the defects studied in Figure |3ji and [3ji 
are quite different, they exhibit a similar trend in the oxidation barriers and in the energy release 
upon oxidation, an encouraging finding that enables a general understanding of the oxidation of 
BP in terms of energy barriers alone, regardless of the specific atomistic arrangements found at 
individual intrinsic defects. 

Geometrical and electronic characterization of the oxidized structures 

We conclude this work with an analysis of the global shape and the electronic properties of the 
oxidized structures that is provided in Figures [ 4 ] and [ 5 } A line defect 30 does not produce significant 
curvature on BP (H ~ 0 and K ~ 0 in Figure[4]), but it creates a periodic compression/elongation of 
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Figure 4: Local shape and electronic properties of the planar BP structure as it absorbs (a) the first 
and (b) the second oxygen dimer. The “structure” subplots represent a segment of the material 
larger than the one seen in Figures |T^,[2ji and [3^. BP has a periodically buckled structure, and the 
geometry is displayed at sublayers .Vj and S 2 along with the local thickness t/to- Oxygen dimers 
(and hydrogen atoms at the finite boundaries of the BP monolayer) gain electrons from BP and the 
electronic gap is reduced as the second dimer is absorbed. Simulated STM images are shown too. 


interatomic distances that is captured by the white/black tones seen on Tr(g) and Det(g) in Figure 
[ 4 } 32 The squares on the individual local shape subplots on this Figure have sides that are 50 A long 
and provide a global view of the shape of the finite systems we work with. 

The chemisorption of a single oxygen dimer does not change the shape of this planar structure 30 
significantly; as indicated in previous paragraphs, BP accommodated the dimer by bringing it out 
of the planar structure and thus maintaining its original shape without any significant alteration (see 
local shape subplots in Figure [4}i). The oxygen atoms receive a large amount of electronic charge 
from phosphorus atoms (see charge transfer subplot in Figure [4ji). Additionally, the hydrogen 
atoms that passivate the edges of the final structure receive electronic charge from phosphorus. 
However, a similar charge transfer from hydrogen atoms is seen throughout all systems studied, 
and for that reason these plots are omitted from now on. 

The next question is whether localized states that decrease the magnitude of the electronic 
gap are created at early stages of oxidation. This question is answered by determining the first 
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electronic wavefunctions below and above the Fermi energy ( Ef ) having a non-zero density at a 
given individual atom. A single chemisorbed CF does not decrease the magnitude of the electronic 
gap (see electronic gap subplot), which remains close to about 1.1 eV in Figure^ due to finite-size 
effects. 32 

Many two-dimensional materials localize electronic states at intrinsic defects. In fact, it is the 
presence of these localized electronic states what permits the facile identification of structural de¬ 
fects. 63 But many structural reconstructions of BP 30 31 lack localized electronic states, making the 
identification of defects under TEM and optical probes a more difficult task. The simulated STM 
images provided at the far right of Figure [4ji are three-dimensional isosurface images produced 
from the density of electronic states having energies up to +1 (down to —1) eV from Ef . 29 64 The 
dislocation line is apparent on the STM at +1 eV. The patterns seen on the STM image at energies 
below Ef are also seen on the electronic gap subplot. All remaining Figures in this manuscript 
were arranged following a layout identical to that of Figure [4ji. 

The electronic properties of the planar BP sample after the chemisorption of a second dimer 
are shown in Figure [4}?. The increase in distances among P atoms once the second dimer was 
absorbed is captured by a new visible dark pattern on the metric invariants in the upper sublayer S 2 
on Figure [4 ]d where the second dimer was absorbed. Curvatures (H and K on Figure [4})) continue 
to be negligible and the ratio t/to remains close to unity, indicating no significant bulging either. 
Importantly, the overall shape observed in Figure [4] tells us that the absorption of oxygen does not 
induce a significant curvature at the onset of the oxidation process. Oxygen atoms continue to 
gain electronic charge as the second dimer is added. A significant difference, though, concerns the 
magnitude of the gap that is reduced this time, as highlighted by the red undertones in the electronic 
gap subplot in Figure [4j>. This localized charge is also evident on the simulated iso-surface plots. 

If the structure has some curvature originally, the oxidation process tends to increase it. To 
demonstrate this statement, the conical structure has a modest protrusion of 4 A and a slant height 
of about 25 A (for the sake of comparison, a graphene cone will have a much larger height equal 
to 14 A for a slant height of 25 A). This variation in height is well within the range of heights 
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Figure 5: Conical structure (a) without oxygen dimers; (b) with one dimer; and (c) with two dimers. 
The shape information now indicates a clear tendency towards amorphization as BP oxidizes: BP 
sees fluctuations on its metric and curvature, and it bulges too as seen on the t/to subplots. Oxygen 
atoms gain electrons from BP, and the electronic gap is reduced as these two dimers are absorbed. 
Simulated STM images are also shown. 


measured for the exposed layer on recently exfoliated black phosphorus stacks. 

The initial shape of the conical structure studied in Figures |TJ), [2jo, and [3j), can be seen on 
Figure [5}i and has the following salient features: the lower (upper) sublayer ,V| ( 52 ) is compressed 
(elongated) near the apex, as indicated by the white (black) features on Tr(g) and Det(g) around 
the structure’s geometrical center. The self-passivated defect is shown in black on Tr(g) and Det(g) 
to the right of the geometrical center, indicating increased metric invariants there. This defect 
did not pin electronic charge within the electronic semiconducting gap, 30 31 as can be seen on the 
electronic gap subplot in Figure [5}l As a result, the band gap remains of the order of 0.89 eV down 
from 1.1 eV on the planar structure in an effect due to curvature. 32 The lack of localized states at 
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the defect 30 31 57 may mask its identification by optical probes. 

The absorption of O 2 in a planar structure did not increase the curvature significantly, but there 
is a clear increase in the local magnitude of the gaussian curvature K in the conical structure, 
Figure [5 }d, as soon as the dimer is absorbed: oxidation increases the curvature of an already non¬ 
ideal structure. In terms of experiments, this means that a structure with a large height distribution 
displays local changes of curvature as well, and the bulges seen experimentally after degradation 
are consistent with an increase of curvature upon oxidation. The gradual increase of curvature is 
reaffirmed in Figure [5]:. The electronic gap of this structure reduces, somewhat dramatically, to 
0.65 eV -signifying a 27% decrease with respect to the gap of the structure seen in Figure [5]i. An 
increase in amorphization 52 is evident when the second dimer is part of the BP structure (Figure 
|5j:), and the electronic properties do not seem to be largely modified from what it was discussed in 
Figure [5J). 

Conclusion 

The search for mechanisms to explain the oxidation of black phosphorus remains one of the most 
exciting and challenging problems in the context of this layered material where chemistry plays 
a fundamental role. We provided a viable mechanism for the photo-oxidation of BP at intrinsic 
defects. Intrinsic defects lower the chemisorption barrier of ideal black phosphorus and make the 
oxidation barrier accessible through the absorption of photons within the visible and ultraviolet 
range, thus enabling the photo-induced oxidation of BP and the dissociation of oxygen dimers. We 
studied energetics of chemisorption and found the oxidation to be highly exothermic. In addition, 
the local geometry, electronic properties, and simulated STM images of these structures were dis¬ 
cussed as well. Further avenues of study concern the degradation of subsequent layers following 
similar mechanisms. This study provides novel insights into the degradation dynamics of this ma¬ 
terial, and gives rise to a plausible mechanism for the oxidation of black phosphorus when exposed 
to light. 
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Methods 


Structural optimizations were carried out with the SIESTA DFT code 6265 66 with spin-polarized 
density functional theory and the PBE exchange-correlation potential, following a conjugate-gradient 
method until all atomic force components were smaller than 0.04 eV/A. All P and O atoms were 
allowed to relax their forces during the optimization process. Movies displaying the chemisorption 
process were created from individual coordinate snapshots as the conjugate-gradient optimization 
went on, following an in-house scripting procedure. The analysis of the discrete geometrical con¬ 
formations is adapted from RefJ 67 (see Ref. 68 too). The “local gaps” are obtained by projecting the 
first state having a non-zero electronic density at any given atom below and above the Fermi level, 
and are indicative of charge localization and level pinning within the electronic gap. STM images 
were created from individual wavefunctions whose density is squared and added up (down) 1 eV 
from T>PH! 

To define the discrete geometry employed on Figures 1, 4 and 5, we consider three directed 
edges ei, e .2 and e 3 such that ej + e 2 + e 3 = 0, and define <2/ = e/ • e/ (/ = 1,2,3), representing the 
square of the smallestyzmYe distance among atoms on the 2-D lattice. 67 



Figure 6: Discrete tensors based on triangulations are expressed in terms of averaged normals n 7 , 
edges e 7 , the normal of the triangle v, and dual edges e* = v x e 7 (/ =1,2,3). 

We consider the change in orientation among normals n j and as well, and project such 
variation onto their common edge e/: This is is, one defines Q 1 / = (n*. — n 7 ) • e/ (see Figure 6; /. k. I 
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are permutations of integers 1, 2, and 3). In previous Equation, ft/ is the average over individual 
normals at triangulated area elements within the polygon surrounding atom l and highlighted by 
dashed lines on Figure 6. The dual edge is defined by ej = v x e/, with v the normal to the triangle 
formed by atoms 1, 2 and 3 and Aj is the triangle area (—ei x e 2 = 2 AjV)W^ 

This way, the discrete metric tensor takes the following form: 32 67 

*=-572 E (Qj-e't- ei)«;®«;, a) 

8A 0 (jXl) 

with Aq the area of the triangulated area element at the reference (non-deformed and defect-free) 
plane. 

The discrete curvature tensor is: 

k = -s 4 E ( 2 " - a" - ef)e* ®e*. ( 2 ) 

6A r (j,k,i) 

The parenthesis (j.k. I ) indicates a sum of three terms, as follows: (j = 1, k = 2, / = 3), (2, 3, 1), 
and (3, 1, 2). Eqns. (??) and (??) become 3x3 matrices with explicit values for Q\, Q 1 / and e* 
from atomic positions. For instance, the discrete curvature tensor has eigenvalues {0, k\ , ki } at 
each triangulated area element, yielding H = (k\ Akj)/! and K = k\ki. The geometrical invariants 
reported at point j are averages over their values at individual triangles sharing this vertex. 

The thickness T is the distance among an atom in the lower sublayer Si and the centroid of 
three nearest atoms belonging to the upper sublayer S 2 . 

Associated content 

Supporting information 

Movies showing the dissociation of oxygen dimers at BP are available free of charge n the ACS 
Publications website. 
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Intrinsic defects, fluctuations of the local shape, and the photo-oxidation of 
black phosphorus 
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Synopsis statement: Intrinsic defects induce fluctuations of a local shape. The 
energy barrier for oxygen dimers to pierce black phosphorus reduces its magni¬ 
tude there within the realm of light-induced excitations. Defects dissociate oxygen 
dimers too. 
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